Cavity Approach to the Random Solid State 
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The cavity approach is used to address the physical properties of random solids in equilibrium. 
Particular attention is paid to the fraction of localized particles and the distribution of localization 
lengths characterizing their thermal motion. This approach is of relevance to a wide class of random 
solids, including rubbery media (formed via the vulcanization of polymer fluids) and chemical gels 
(formed by the random covalent bonding of fluids of atoms or small molecules). The cavity approach 
confirms results that have been obtained previously via replica mean-field theory, doing so in a way 
that sheds new light on their physical origin. 

PACS numbers: 



Introduction — Permanent random chemical bonds, when 
introduced in sufficient number between the (atomic, 
molecular or macromolecular) constituents of a fluid, 
cause a phase transition — the vulcanization transition — 
to a new equilibrium state of matter: the random solid 
state. In this state, some fraction of the molecules are 
spontaneously localized, and thus undergo thermal fluc- 
tuations about mean positions, the collection of which 
positions are random (i.e. exhibit no long-range regular- 
ity). As this is a problem of statistical mechanics (the 
constituents are undergoing thermal motion) in the pres- 
ence of quenched randomness (the constraints imposed 
by the random bonding), it has been addressed via the 
replica technique. Indeed, the example of cross-linked 
macromolecular matter served as early motivation for 
Edwards' development of the replica technique (see, e.g., 
Ref. 0). 

Our purpose in this Letter is to consider two central 
diagnostics of the random solid state: the fraction Q of 
localized particles and the statistical distribution A/"(£ 2 ) 
of squared localization lengths £ 2 of these localized parti- 
cles. Specifically, we show how results for these quantities 
can be obtained in an elementary and physically trans- 
parent way, via the cavity method. The cavity method 
has proven flexible and powerful in the analysis of a va- 
riety of other disordered systems, e.g., spin glasses 0. 
The present work is based on the version used to address 
spin glasses having finite connectivity 0. 

The results that we shall obtain via the cavity method 
are amongst those already known via a (less elementary 
and less physically transparent) application of the replica 
technique, together with a mean-field approximation; for 
reviews see Refs. 00. Thus, it was already known 
that (as is typical for mean-field theories) Q obeys a tran- 
scendental equation, in this case 



1 -Q = exp(-^ 2 Q), 



(1) 



where /i 2 is a parameter that controls the density of 
cross-links. The instability of the fluid state and its re- 
placement by the random solid state are signaled by the 
emergence, as \x is increased beyond a critical value (here 



unity), of a positive solution to Eq. QJ, although the 
formulation is not restricted to the critical regime. This 
result for Q was, in essence, found by Erdos and Renyi, 
in their classic work on the statistical properties of ran- 
dom graphs pfl. As for the distribution A/"(£ 2 ), at the 
mean-field level and for near-critical values of the cross- 
link density, the replica approach yields the scaling form 



(2) 



where e [oc (/i 2 — 1)] measures the distance from the crit- 
ical point, and in which the scaling function tt obeys the 
nonlinear integro-differential equation 



de b 



ir(9)) =tt(0) - (vtott) (6>), 



(3) 



where tt o tt indicates a Laplace convolution [9j . This re- 
sult was obtained from a semi-microscopic approach by 
Castillo et al. and re-derived via a Landau-type the- 
ory by Peng et al. ^(j- Some support for the results for 
Q and tt has been obtained numerically pl| and experi- 
mentally [T^ . 

More recently, the mean-field level results have been 
improved in two directions: via renormalization-group 
analysis in the vicinity of the upper critical dimen- 
sion (which is six for this random solidification transi- 
tion) [T10, and via analysis of the Goldstonc fluctua- 
tions | 15| (notably in two dimensions, where their conse- 
quences are — not surprisingly — dramatic). We pause to 
mention that we view the cavity method as complement- 
ing rather than supplanting the replica method: the lat- 
ter provides access to the powerful array of field-theoretic 
tools, the former opens the way to a perturbative treat- 
ment of correlations beyond mean-field. 

Cavity method for randomly cross-linked macromolecules, 
especially in the vicinity of the random solidification tran- 
sition — We begin by considering a system of vulcanized 
macromolecules, as depicted in Fig. QJi. We characterize 
the system by the fraction Q of localized chain segments 
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FIG. 1: (a) Snapshot of a cross-linked (open circles) macro- 
molecular system, in which a fraction Q of the macromolecules 
belong to the infinite cluster (and are thus localized), and the 
remaining fraction 1 — Q are delocalized. A localized segment 
exhibits thermal position-fluctuations over a length-scale £ 
(i.e. its localization length). Due to the heterogeneity of the 
network, distinct segments have distinct localization lengths, 
(b) Coupling of a new macromolecule (dashed) to the system 
in order to probe the distribution 7V(£ 2 ). The new macro- 
molecule has N c contact points with the existing network, 
and cross-links are introduced independently at these points. 



and the statistical distribution M of mean squared lo- 
calization lengths £ 2 of the localized particles. We then 
envisage adding a further macromolecule to the system, 
as shown in Fig.njj. Of all the segments on this chain, we 
suppose that a certain number N c are sufficiently close 
to segments of the original system to have a chance of 
becoming cross-linked to them. We suppose that fluctu- 
ations in this number are sufficiently small that we may 
neglect them. Next, we consider a random cross- linking 
process that results, with independent probability p, in 
cross-links actually being introduced between each of the 
N c close pairs. Within this framework, the probability 
that exactly k cross-links are introduced is then given by 
the binomial formula: ( N k ")p k (l — p) Nc ~ k . We now ask 
the question: What is the probability Vk that exactly k 
of these cross-links are made to localized segments? To 
answer it, we make the approximation that the probabili- 
ties for the segments of the original system to be localized 
describe independent random variables, in which case the 
probability of any one such segment being localized is Q. 
Then, collecting together the contributions to this proba- 



bility, which arise from k' (with k < k! < N c ) cross-links 
being formed, of which exactly k are to localized seg- 
ments, we arrive at the formula 

^ = E ( N ;)p k \^p) N - k '( k ')Q^^Qf- k - (4) 
k'=k \ ' ^ ' 

Via a straightforward application of the binomial theo- 
rem one can perform this summation, and hence arrive 
at the result 

Vk=( N k c y P Q) k (l- P Q)^- k . (5) 

Let us evaluate these probabilities for the three cases of 
relevance, viz., 

V = (1- P Q) N % (6a) 
V\ = N cP Q{l-pQ) N *-\ (6b) 
V 2 = \N c {N c -l){pQ) 2 {l-pQ) N <- 2 . (6c) 

In the limit of main interest, viz. pQ < 1, these proba- 
bilities simplify to 

(r 0) Vi,V 2 ) « e- N ' pQ (l,N cP Q, ^N C (N C - \){ P Q) 2 ). 

(7) 

To arrive at a self-consistent equation for Q (as a func- 
tion of p and N c ) we require that the probability of the 
added macromolecule being cross-linked to exactly zero 
localized segments be 1 — Q, which gives 

l-Q = (l-pQ) N *. (8) 

In the limit p Q <C 1 this becomes 

l-Q = cxp(-^ cP Q), (9) 

i.e., Eq. iJTJ, provided we make the (physically sensible) 
identification p N c = n 2 . 

To arrive at a self-consistent equation for Af (as a func- 
tion of p, N c and £ 2 ) we address the probability that the 
segment at arc-length a (with < a < 1) on the added 
macromolecule has squared localization length £ 2 . As the 
segments of the added macromolecule will only be local- 
ized if they are attached to at least one localized segment 
of the original system, we should replace the probabilities 
{■Pfe} 2 =0 by the probabilities conditioned on at least one 
attachment being to a localized segment of the original 
system (i.e. k = 1,2). Thus we arrive at the probabilities 
(V\, V-i) which we write in the form (1 — a, a). Consider 
the situation in which the added macromolecule is at- 
tached at its arc-length g\ to a single segment of the 
original system, that segment having squared localiza- 
tion length £ 2 , as depicted in Fig.^t- Furthermore, sup- 
pose that the added chain is Gaussian. Then, by the el- 
ementary properties of random walks, the mean squared 
spatial separation of segments separated by arc-length 
|<7 — <7i| is given by |er — a% \ 1Z 2 where 1Z 2 is the mean 
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FIG. 2: (a) The added macro-molecule is attached at a sin- 
gle segment to the infinite cluster (i.e. the shaded region), 
(b) The added macromolecule is attached at two segments to 
the infinite cluster. 

squared end-to-end distance of each chain. Thus, if Af is 
the distribution for the squared localization length of the 
segment of the original system to which the new chain is 
attached at arc-length o\ , the distribution of the squared 
localization lengths £ 2 for the segment at arc-length a on 
the new chain will be given by 

J de 1 m 2 i)s(e-(e 1 + w-<yi\ni)). ao) 

Now, supposing that the addition to the squared local- 
ization length, | (7 — <7x|"7£g, is small compared with the 
localization lengths that feature with appreciable weight 
in Af, this approximates to 

n (e -\a- oii^) « m 2 ) -w- <Ti\niu'{e)- 

Next consider the situation in which the added macro- 
molecule is attached at its arc-lengths o\ and a 2 to two 
segments of the original system, these segments having 
respective squared localization length £ 2 and £f , as shown 
in Fig. [2b- In fact, as the probability that the added chain 
has two cross-links to localized segments of the original 
system is (in the limit of interest) much smaller than the 
probability that it has one cross-link, we shall not need to 
keep track of the arc-length locations of the cross-links; 
in this situation it will be adequate to treat the added 
chain as a point object. Then, as this object is attached 
to two localized objects, it is pinned more sharply than ei- 
ther, this parallel form of pinning giving rise to a smaller 



squared localization length £ 2 , via the formula 

r 2 = er 2 + ^ 2 - (12) 

So, assuming that £ 2 and £f are independent [and thus 
governed by the joint distribution Af(£,f) A/"(£f)L the dis- 
tribution of £ 2 is given by 

/ de 1 mi)de 2 me 2 )s(e-(^ 2 +^ 2 y 1 )- w 

We now put these results together to construct the 
distribution of squared localization lengths for segment 
a of the added chain, arriving at 

{i-a)(N(e)-\o-vi\KN'{e)) 
+ a Jde 1 mi) dt,imt) s(e (er 2 + s-y 1 ) (14) 

Finally, we average over the segment a of the added 
chain, as well as the location o\ of the cross-link, us- 
ing J da da' \a — a 1 | = 1/3, thus arriving at the self- 
consistent equation obeyed by the distribution of squared 
localization lengths: 

m 2 ) = a-a)(m 2 )- hK^'tt 2 )) 

Observe that by integrating both sides over £ 2 and in- 
voking the property that Af vanishes at the limits £ 2 = 
and oo, we recover the normalization condition that 
/ g?£ 2 A/"(£ 2 ) = 1. The scaling property of Af shows up 
via the following change of dependent and independent 
variables: 



2 1 — a /ci 

under which Eq. (|15(l becomes the sought integro- 
differential equation J2J). We identify the parameter e in 
Eq. J2J with rrj^R • Notice that the cavity method al- 
lows us to compute corrections to Eq. @ perturbatively 
in 1/Hg. 

Cavity method for randomly bonded Brownian particles 
at arbitrary bonding densities — The cavity approach can 
be extended to address the chemical gelation transition, 
i.e., the transition triggered by the introduction of ran- 
dom covalent bonds between atoms or small molecules 
(rather than macromolcculcs) in the liquid state. Wc 
shall address the model studied previously by Broderix 
et al. 0| , which consists of a collection of point particles 
undergoing Brownian motion at a certain temperature. 
Permanent bonds are then introduced at random between 
nearby particles, so that pairs of bonded particles be- 
come constrained softly (i.e. by a spring-like, harmonic 
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potential), so that the probability distribution <j> of their 
separations — r k is Gaussian and characterized by a 
length-scale I; 



(r,- - v k ) oc cxp (-|rj - Y k \ 2 /2l 2 ) 



(17) 



where r^/. are the position vectors of particles j and k. 
This model is sufficiently simple that the analysis of it 
need not be restricted to the critical region. 

To approach the statistics of this system of randomly 
bonded Brownian particles using the cavity method, we 
consider the process of adding a new particle. The com- 
binatorics of the bonding follows the form taken for the 
system of cross-linked macromolecules; we simply need 
to convert the notion of contact points into a spherical 
region of a certain radius in which the likelihood of par- 
ticles being bonded to one another is concentrated. This 
sphere is centered on the new particle and, on average, in- 
cludes -ZV C of the existing particles. Then, bonds are ran- 
domly introduced, with probability p, between the new 
particle and some of the N c existing particles that are 
nearby. Thus, the foregoing combinatorics continues to 
apply, and we arrive at the formula for the probability of 
having exactly k bonds with the infinite cluster given in 
Eq. (JSJ. As a consequence, we obtain the foregoing result 
for the fraction of the infinite cluster Q, Eq. 

The physics of the localization lengths is, in fact, sim- 
pler for bonded Brownian particles. When the new par- 
ticle is connected via a spring of length-scale I to one 
localized particle, its localization length Si is given by 
E^ 2 = (£ 2 + i 2 ) -1 . When it is connected in parallel via 
identical springs to k localized particles, its localization 
length Sfc is given by 



(18) 



^k 



tion length of the new particle, we shall average over 
all possible numbers of bonds, weighted with their cor- 
responding probabilities. These probabilities follow from 
the probabilities Vk given by Eq. JSJ, but normalized by a 
factor Q _1 because the new particle will only be localized 
if it is bonded to at least one particle in the infinite clus- 
ter. Hence, we arrive at a self-consistency equation for 
localization-length distribution for the randomly bonded 
Brownian particle model: 



fc=i j ° 



as) 



To construct the distribution of the squared localiza- 



where the conditional probabilities Vk are given by Vk — 
Vk/Q (for k — 1, 2, 3, . . .). The distribution of localiza- 
tion lengths for the randomly bonded Brownian particle 
model was studied previously by Broderix et al. [l|| using 
the replica method and a Mayer cluster expansion. To see 
that the cavity approach result (|19f) recovers their result, 
we take the limit iV c — > oo whilst keeping finite the mean 
number of bonds from the new particle (either to the in- 
finite cluster or to delocalized particles), i.e. pN c . In this 
limit, the binomial distribution tends to a Poisson distri- 
bution: V k — >(0pN c ) k Q k - 1 /W)e- pN c Q . The result of 
Broderix et al. [l|| then follows from Eq. I|19|) by (i) trans- 
forming to a distribution for r = l/£ 2 , and (ii) making 
the identifications k = l~ 2 and c = pN c (i.e. the mean of 
number of bonds associated with a single particle). 
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